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Using continuity, we derive a renormalized Hamiltonian from the parent t-J model to describe 
the properties of underdoped cuprates. The theory is constrained to agree with the behavior at half 
filling, which is well described by the Arovas-Auerbach valence-bond state in which bosonic spinous 
are paired into singlets. Spinon states evolve continuously into the doped region preserving their 
symmetry. We assume that moving holes rapidly destroy magnetic order, which leads to a gap in the 
spinon spectrum and strongly renormalizes the theory. The spin gap leads to two new types hopping 
terms for renormalized holes. In one, a fermionic holon hops within the same sublattice accompanied 
by a singlet backflow, giving rise to a non-Fermi liquid normal state with novel properties. Spinon 
singlets condense below a pseudogap temperature T* (< spin gap temperature r°), which allows 
holons to propagate coherently, forming a spinless Fermi liquid, but without an observable holon 
Fermi surface. Above T*, holons are localized. This is the so-called strange metal phase, which 
is actually a new type of insulator. In the second term a pair of holons belonging to opposite 
sublattices hop, accompanied by a singlet backflow. In the presence of the singlet condensate holon 
pairs condense, leading to d-wave superconductivity; the symmetry is primarily determined by the 
symmetry of the valanece-bond state at half filling. The metal and the superconductor preserve 
the two-sublattice character of the valence-bond state. A careful examination of the nmr, tunneling 
and transport data shows that the predictions of the theory is consistent with experimental results. 
Remarkably, the existence of the spin gap provides a natural explanation for the phenomenon of two- 
dimensionality of the normal state in the presence of interplane hopping. The marked asymmetry 
between hole-doped and electron-doped cuprates is also easily explained. 
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I. Introduction 

The origin of high-temperature superconductivity in 
doped cuprates [l| is believed to be closely linked with 
the unusual behavior of the normal state, which is not a 
Fermi liquid. Furthermore, in the underdoped region a 
pseudogap appears below a temperature T* , much above 
the superconducting temperature Tc Even before 

the discovery of the pseudogap phase, Anderson argued 
that the non-Fermi liquid behavior is due to the two- 
dimensional (2d) nature of the normal state @, and its 
proximity to the undoped phase, which is a Mott insu- 
lator, or equivalently, a quantum antiferromagnet. The 
2d behavior is unexpected, and its origin is not under- 
stood. Anderson also proposed the resonating-valcncc- 
bond (RVB) state, in which spins are paired into singlets 
within the insulator 0]. Upon doping with holes, this 
state would evolve continuously into the doped region 
where spinless holons, which carry charge, would propa- 
gate coherently and create a metal. The elementary ex- 
citations are holons and spin- 1/2 spinous. The physical 
electron is a composite particle. The expectation is that 
the apparent spin-charge separation would lead to non- 
Fermi liquid behavior, and further that the spin singlets 
would acquire charge via some interaction with holons 
and become superconducting. 

Theoretical studies are usually based on the large- [/ 
Hubbard model, or, equivalently, the t-J model on a 
square lattice, where J = /U is the exchange coupling 
between spins, and t describes nearest-neighbor hopping 
of electrons such that no site is doubly occupied. For 



cuprates, t/J ^ 3 — 4. The constraint of no double occu- 
pancy corresponds to a U{1) gauge symmetry Q. Much 
of the work has been devoted to deriving an effective 
low-energy theory of propagating spinons and holons (or 
pairs), coupled to a gauge field. In principle, such a task 
should not be too diffcult since one can invoke continuity. 
The effective action would depend on the symmetries of 
the underlying vacuum state, i.e., a renormalized version 
of the RVB state of paired spinons, which is presumed 
to describe the insulator at half filling. After more than 
two decades of intensive work, a successful theory has not 
been found. 

This failure is surprising since, using a mean-field ap- 
proximation, Arovas and Aucrbach [5| have shown that 
the behavior at half filling is indeed reasonably well de- 
scribed by an RVB state in which singlets formed by pair- 
ing bosonic spinons (Schwinger bosons) condense below 
a temperature Tuvb- (In this paper we use the term 
RVB to describe any valence-bond state) . One complica- 
tion is that a fraction of spinons remains unpaired, and 
condenses separately to give rise to long-range antiferro- 
magnetic order. This is not a defect of the theory since 
it correctly describes the ground state which is known 
to be ordered @. In other words, the system is in a 
mixed phase of an RVB state and a Nccl state, with 
Trvb > Taf (hi d = 2, Taf = 0). In cuprates the AF 
insulating state exists up to a hole density x of about .05, 
beyond which there is a transition to a pseudogap metal 
(and a d-wave superconductor), with no long-range mag- 
netic order of any type. However, a seemingly straightfor- 
ward extcntion of the Arovas-Auerbach MF theory (with 
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fermionic holons) to the doped region leads to a metallic 
state with spiral magnetic order [7[ . This failure has lead 
to a virtual abandonment of the approach based on the 
Arovas- Auerbach RVB state and bosonic spinons Q . 

Instead, RVB theories based on the slave-boson repre- 
sentation (spinons are fermions, holons are bosons) have 
been used widely [§, [13, [U (for a review, see Lee et 
al. But these do not work at half-filling, and to 

our knowledge, there are no experimental signatures of a 
transition to a different RVB state accompanied by sta- 
tistical transmutation upon doping. There are serious 
doubts as to their efficacy in describing the physics of 
underdoped cuprates. Indeed the failure to connect with 
the physics at half-filling is one of the central problems 
in high- Tc theory. 

We have solved this problem by deriving a renormal- 
ized theory for small x which is consistent with the 
physics at half filling, but does not have the problems 
of spiral instability. The key point is that, for the t-J 
model, the spiral state has been shown to be unstable 
p^ . We therefore assume that (1) the spin states evolve 
continuously from half filling, preserving their symmetry 
and particle statistics, and (2) moving holes rapidly de- 
stroy long-range magnetic order at T = beyond some 
small critical concentration Xc, strongly renormalizing 
the theory in the process; but the VB state survives up to 
a temperature T]iyB{x). These are reasonable assump- 
tions, and formed the basis of early gauge theories based 
on the Schwinger-boson representation [l3 . [ist . But in 
those theories it was assumed that the initial renormal- 
ization leads to a sublattice-preserving t'-J model, where 
t' is an effective next-nearest neighbor hopping parame- 
ter. 

The actual renormalized Hamiltonian is found to be 
quite different. A brief account of the theory focusing 
primarily on the pseudogap phase and some of the re- 
sults have been reported earlier (l6j . Here we provide 
a detailed description, and show that it not only ac- 
counts for the main features of the phase diagram, but 
predicts novel properties of the underlying states that 
require a major revision of our current understanding 
of the cuprate phenomenology. Additionally, the theory 
provides a natural explanation for the two-dimensionality 
of the normal state, as well as for the marked asymmetry 
between hole-doped and electron-doped cuprates. 

The continuity requirement ensures that there are ex- 
actly the same three spinon states in the doped region 
as at half filling, and holon motion do not introduce any 
new order parameters for the spinons. Destruction of 
magnetic order automatically leads to a gap for bosonic 
spinons in the other two states: the RVB state, and the 
nonordered state [1, ■ The spin gap - not to be con- 
fused with the pseudogap - is the distinctive feature in 
the theory. As described in section IV, the gap allows us 
to renormalize away nearest- neighbor hopping (i), and 
use perturbation theory to obtain a minimal Hamiltonian 
involving sublattice-preserving hopping of renormalized 
holes, accompanied by singlet backflows. Two types of 



new hopping terms are generated: in one a single holon 
hops within the same sublattice, and in the other a pair 
of holons hop together. 

We analyze the Hamiltonian in section V, showing 
that the one-hole term gives rise to exactly two 'normal' 
states, as found in cuprates, whose symmetries are de- 
termined by the underlying spinon states (at half filling) 
- no new order parameter is introduced, either in the 
spin or the charge sector. In the absence of singlet con- 
densation (i.e., above Ti^vb{x)); holons are localized. We 
identify this state with the so-called strange metal phase, 
except that it is not a metal, but a new kind of insula- 
tor, which is an important prediction. Coupling to the 
RVB condensate allows holons to hop coherently within 
the same sublattice, and create a spinless Fermi liquid of 
concentration x below TuvBix), which we identify with 
the pseudogap temperature T*. However, gauge invari- 
ance ensures that there is no observable (small) Fermi 
surface, which is another prediction. 

As shown in section VI, in the presence of a RVB con- 
densate the pair-hopping term allows holon pairs to Bose 
condense, giving rise to d-wave superconductivity. The 
symmetry is largely determined by that of the underly- 
ing RVB state known from half filling. As discussed in 
section II, the valence-bond state has a two-sublattice 
property of its own: on average, singlets connect spins 
on opposite sublattices. This property is predicted to be 
preserved in the pseudogap metal and superconducting 
phases. 

The absence of any new order (other than supercon- 
ductivity) imposes constraints on the number and type 
of phases; the resulting phase diagram is in qualitative 
agreement with the experimental one. The constraints al- 
low us to extract many universal features from the struc- 
ture of the Hamiltonians in various phases, without de- 
tailed calculations or fine tuning. For example, the spin 
gap by itself rules out a conventional Fermi liquid state 
for small x. In section VII we carefully examine the ex- 
isting experimental results, and find strong support for 
the predictions of our theory. In particular, results from 
nmr, transport, optical conductivity and tunneling mea- 
surements, taken together, provide strong evidence for 
a spin gap, a nonmetallic phase above T*, and a spin- 
less Fermi liquid without a Fermi surface below. We will 
point out that these results can not be reconciled with 
either a Fermi liquid, or earlier RVB theories. 

In section VIII we discuss the fundamental issue of 
the two-dimensionality of the normal state, which is un- 
expected since the cuprates are 3d - though layered - 
materials, and localization within a plane would cost too 
much kinetic energy. Indeed, all the other nearby states: 
the antiferromagnet in the undoped phase, the supercon- 
ducting state, and the Fermi liquid state that appear at 
large doping exhibit 3d behavior. Hence theories that 
are based on the 2d model are suspect if they fail to ex- 
hibit 2d confinement when hopping in the perpedicular 
direction is turned on. Indeed, such a test will probably 
invalidate most of the existing theories. We show that, in 
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our case, the spin gap provides a protective mechanism 
against delocahzation of holons in the z direction, which 
is not obvious since holons arc spinless. However, the 
holon pairs can hop coherently, which is essential for the 
observed 3d superconductivity. We solve the 3d problem 
and show that condensation energy is enhanced due to 
3d coupling. 

Another issue that is not completely understood is 
the significant difference between the hole-doped and 
electron-doped cuprates. In section IX we will argue that 
the difference is due to direct intra-sublattice hopping 
t', which breaks electron- hole symmetry in the Hubbard 
model. Normally the difference is not significant, we will 
present results to show that in our case it is substan- 
tially enhanced becuse the nearest-neighbor hopping (t) 
is renormalized away. The issue is subtle and depends on 
the symmetry of the underlying spinon states. 

Our renormalized Hamiltonian belongs to the class 
of short-range RVB models considered by Kivelson, 
Rokhsar and Sethna [H, [11]. However, because of 
the consistency requirements the Hamiltonian differs 
significantly from those proposed earlier. The present 
theory is based on a paper published much earlier, 
in which a self-consistent perturbation expansion in 
powers of hopping was used to study the destruction of 
long-range magnetic order and the occurrence of non- 
Fermi liquid behavior [20| . The intra-plane pair hopping 
mechanism for superconductivity in the present theory 
was also discovered there. An effective Hamiltonian for 
renormalized holons was derived in ref. [2l] |. but with 
an (assumed) wrong (spiral) normal state. In section 
ni we review these and earlier works by other authors 
that are relevant to the present derivation. We draw 
attention to the early single-hole results which show 
clear evidence for strong renormalization, and effective 
sublattice-preserving hole hopping. Our conclusions are 
summarized in section X. 

II. Model and Symmetries 

The physics of no double occupancy is taken into ac- 
count by representing the electron as a composite object 
created by = hi , where fej^. creates a spinon of spin 
a and hi destroys a holon at lattice site i, subject to the 
constraints that number of holons plus spionon at each 
site is one. Spin operators are represented as = 6^|6ix, 

and SI = ^(&||&it ~ ^ij^ii)- Then the t-J Hamiltonian is 
given by 

U^-tY^clc,. - 2J^Al,A,„ (1) 

where the sum is over nearest neighbors. Ay — \ — 
6i|6j|] destroys a singlet connecting spinous on nearest 
neighbor sites. The second term describes exchange in- 
teraction written in terms of Ay's, which makes the role 
of the singlets explicit. For the most part we consider 
the square lattice. However, later we will need to include 



next-nearest-neighbor hopping within the plane and as 
well as hopping between planes. 

The Hamiltonian is invariant under a local UiX) gauge 
transformation [i^]: 

^IG ^ ^ifT^ ; hi(j > hifjC . (2) 

The number of spinous plus holons 

M = bl^bi^ + bl^bii + h\h,. 

is then conserved at each site. The physical model lives 
in the projected subspace defined by A/i = 1, for each i. 
We can choose any statistics for spinous and holons as 
long as the statistics of the gauge invariant quantities are 
correctly given. However, in an effective theory, in which 
spinous (and holons) are the low-energy excitations, we 
expect the system to choose a particular statistics. 

Half-Filled Case: 

At half-filling the hopping term can be dropped. In the 
mean-field approximation constraints are treated on the 
average, and singlets condense so that the valence-bond 
order parameter Aij =< Ay > becomes nonzero. For 
dimensionality > 2, the MF theory works well at low 
T if spinous are bosons. For the symmetry of the order 
parameter we can choose: 

A,^. ==yle'^'5-('-'-'-j), (3) 

where Q is the zone corner wave vector. In d = 2. Q = 
(tt, tt). However, for c? > 2 it has been shown that the 
ground state is not a pure RVB state since a fraction of 
spinous remains unpaired 01 • Since spinons are bosons 
they condense independently so that < bi^ >^ 0. This 
leads to a nonzero value for < Sf >=< bj-^bi^ >, giving 
rise to a two-sublattice antiferromagnetic order (in the 
x-direction). 

The ground state is therefore a mixture two competing 
phases, characterized by distinct order parameters, in 
which magnetization is reduced from its classical value 
by the presence of the singlet condensate. The MF 
theory describes the ground state properties in d = 2 
rather accurately, and in o? = 3, reproduces the results of 
Holstein-Primakoff quantum spin- wave theory Q . This 
success is remarkable, but not accidental, since the MF 
equations are identical to the equations describing the 
spin- wave theory In fact, the latter can be derived 
from the Schwinger-boson theory by integrating out 
the constraints, and expanding the action in powers of 
1/5. At finite T, AF order disappears above Taf (hi 
d = 2, Taf = 0), but the RVB state survives up to 
Trvb > Taf- 

Symmetry of the RVB State: 

In the MF approximation, condensation of singlets 
(and spinons) breaks gauge symmetry. Similarly, mo- 
bile holons lead to < hjhj >^ for i ^ j, which also 
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breaks gauge symmetry. However, such local symmetries 
can not be spontaneously broken since there is no 
rigidity with respect to phases of the order parameter 
Aij generated by the transformation (2). When aver- 
aged over these gauge equivalent choices, Aij (and other 
gauge- variant quantities) vanish. In this sense there is no 
condensation. Similarly, spinon and holon Green's func- 
tion are not gauge invariant, so these particles are not 
directly observable. The remaining phase fluctuations of 
the order parameter are studied via a gauge theory, as 
has been done by many authors. However, there is still a 
transition across Trvb since the free energy is gauge in- 
variant, and retains the characteristics of the condensate. 
In the uncondensed state (i.e., above Trvb) both phase 
and amplitude fluctuations are strong, and spinous and 
holons remain localized. In the condensed phase, fluc- 
tuations are weak. In this paper we will use the term 
condensation in this sense, and the term symmetry of a 
MF state to reflect the symmetry of the whole gauge- 
equivalent class. 

Furthermore, condensation can lead to certain type of 
gauge-invariant order on the lattice. For example, the 
RVB state has an important two-sublattice property (see 
below). Other choices of symmetry lead to flux phases. 
Similarly condensation of spinon pairs and holon pairs 
lead to condensation of gauge-invariant composite enti- 
ties such as a electron pair, whose symmetry depends on 
the symmetry of the constiuents. 

Although gauge fluctuations arc supposed to be weak 
in the ordered phase at low T, for a detailed quantita- 
tive theory (particularly for transport properties) one has 
to include them [2^. Similarly, the MF theory will not 
work well near Trvb where it will be necessary to con- 
sider both amplitude and phase fluctuations. For the 2d 
system these may be strong enough to convert the appar- 
ent second-order MF transition to a Kosterlitz-Thouless 
type transition (see [HI). We will not study these issues 
here. 

At half filling, the MF theory has been formally ex- 
tended by Read and Sachdev to situations where AF or- 
der is destroyed at T = by quantum fluctuations pTf . 
In the ordered phase spinous are gapless, leading to gap- 
less spin-wave excitations. However, once magnetic order 
is destroyed, a gap A., appears in the spinon spectrum. 
Both cases can be described by the same MF theory with 
different values of A, which acts as the spinon chemical 
potential. The spinon spectrum has the form 

c.(k) = [A2-0(k)2]^ (4) 

where 

(/)(k) = 4 J^(sin k,j. + sin ky) (5) 

is the spinon 'gap' function. The minima of tj(k) arc at 
±(7r/2, 7r/2). In the ordered state, A is chosen so that 
(x'(k) vanishes at these points leading two-sublattice AF 
order. 

In the absence of AF order, spinous are gapped since 
^min > 0. Then, as shown by Read and Sachdev gauge 



fluctuations lead to a proliferation of instantons, \yhich 
results in the confinement of spinous into pairs [l7j . In 
this case the Arovas-Auerbach type state is unstable 
relative to a valence-bond solid. In our case, we do 
not need to worry about this, since there is no gap due 
to the presence of AF order. Away from half filling 
a spinon gap does appear, but, as we shall sec, holon 
motion causes singlets themselves to hop around, which 
would destabilize the VB solid. 

A. Correlation Functions: 

Like the Neel state, the bosonic RVB state has a two- 
sublattice character of its own. From Eq.(5) we see 
that although 0(k) is not gauge invariant, it satisfies the 
gauge-invariant condition: 

0(k)=0(Q-k). (6) 

Now consider the spinon correlation function Bjm = 
i < h\^^bju > for any two sites j and m. In the 
MF approximation this is given by 

= ^Ec°«(l^-('^"-^j))^[l/2 + ?^(^(k))], (7) 

where n(w) is the Bosc function. Using Eq.(5), we find 

= Bj„i cos(Q.(rm - rj)). (8) 

Therefore, Bjm vanishes if j, m are on opposite sublat- 
tices. For j,m on the same sublattice Bjm is nonzero if 
the singlet condensate exists {A ^ 0). We will call this 
behavior even. Similarly, consider the anomalous corre- 
lation function, Ajm = \ < (^jT^mi ~ ^iJ-^mt) >i again 
defined for any two sites j, to. (For nearest neighbors this 
is the RVB order parameter). It is given by 




From which we find 

Ajm = -Ajm cos(Q.(rin - rj)). (10) 

Thus Ajm is zero if j, to are on the same sublattice. For 
J, TO on opposite sublattices, Ajm is nonzero as long as 
the singlet condensate exists. We will call this behav- 
ior odd. These symmetry properties arc intrinsic to the 
RVB state since they hold even after long-range AF or- 
der is destroyed, and for any T < Tf/vB- They are also 
gauge invariant, even though the correlation functions 
themselves are not, and therefore have observable conse- 
quences. Thus the spin-spin correlation function is given 
by Ssp.tj =< S+SJ- >= + which, as ex- 

pected, alternates in sign. 

These properties imply that on average RVB singlets 
connect spinous on opposite sublattices. Therefore, 
the bosonic RVB state is similar to the short-range 
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RVB state considered by Kivelson, Rokhsar and Sethna 
[T§ |. What we have shown here is that this RVB state 
has the same red-bhic property even in the presence 
of long-range magnetic order, i.e, even when magnetic 
correlations are not short-ranged. 

B. Magnetically Disordered Phase 

Away from half-filling quantum fluctuations due to 
hole hopping is expected to destroy long-range AF or- 
der at T = 0, leaving behind the RVB phase. Then 
bosonic spinous acquire a gap A^, which is related to the 
magnetic correlation length • Near its minima the 

spinon energy has the form: a;(k) w [A^ -I- Cgfc^], where 
Cs = 2JA '-^ TuvB, and k is measured relative to the min- 
ima. Then the correlation functions Bij and Aij behave 
as 

e-'-/2?cos(iQ.(rj-ri)), 

A,,^-e--/2fsin(iQ.(rj-rO), 

where = jr^ — rj|, and ^ = 2cs/As is the spin-spin 
correlation function. 

The spinon gap, and the consequent exponential 
decay of AF correlations is the distinctive feature of the 
bosonic RVB state. In contrast, spinous are gapless in 
a femionic (i.e., slave boson) RVB state, which would 
lead to power-law correlations. Similarly, a Fermi liquid 
also has gapless spin excitations. Also, in the latter 
two cases, the magnetic correlations are peaked at a 
concentration-dependent wave vector q, which equals Q 
only at half filling. 

III. Aviray from Half Filling: Early Work 

A. One-hole Physics and the t' — J Model 

Several authors have studied the behavior of a single 
hole moving in a magnetic background by expanding in 
powers of hopping, and summing a selected a class of 
diagrams [l^, Ull- R has been found that because of 
spin mismatch the hole can not hop coherently on the 
opposite sublattice, but it can move coherently within 
the same sublattice. For large t/J, the bandwith corre- 
sponding to the coherent peak was found to scale with 
J, reflecting strong renormalization of holes by spin fluc- 
tuations which makes the renormalized hole large and 
heavy. The short-range incoherent hops due to t leads 
to an incoherent spectral background of order t. These 
results are supported by exact numerical work on finite 
lattices containing one hole [2^ . They can be used to es- 
timate renormalized parameters for a many-hole theory. 
However it is erroneous to construct a low-energy theory 
based on the one-hole spectrum, since a single hole does 
not destroy magnetic order the presence of which leads 
to a pole in the electron's Green function. 



Motivated by the single-hole results, several authors 
have developed gauge theories for the many-hole system 
[13, [ill based on the Schwinger boson RVB states (thus 
hoping to establish a connection with the Mott phase), 
with the assumption that the short-range incoherent 
hops caused by inter-sublattice hopping {t) destroys AF 
order and strongly renormalizes the theory, leading to 
a sublattice-preserving t' ~ J model, where ~ J is a 
renormalized next-nearest- neig hbor hopping parameter 
for the renormalized holes [1^. While the assumptions 
are correct, the actual renormalized Hamlitonian, as we 
show later, is quite different. 

B. Spiral and Other Nearest-Neighbor States 

One reason why Schwinger boson methods have not 
been taken seriously is the fear of a spiral instability. 
For a; > 0, moving holes are expected to rapidly destroy 
long-range AF order. As a physical hole hops, it car- 
ries its spin, creating a string of wrong (ferromagnetic) 
bonds, which cost energy and can not be easily repaired. 
In cuprates the AF insulating phase disappears beyond 
a small Xc, giving way to a metal with no long-range 
magnetic order. Now, we can obtain a metallic state by 
extending the mean-field approximation to the hopping 
term [3| , which can be written as 

-t^[BytDij+/i.C.], 

ij 

where Dy = h|hi and By = ^ E<t bj^bia- In the MF 
approximation, in addition to singlets and spinons, the 
composite bosons created by By and Dy also condense, 
so that the hopping term becomes 

ij ij,a 

where B*j =< By > and Dij Dy >, are the av- 
erage spinon and holon backfiows, respectively. Appro- 
priate choices of symmetry for these order parameters 
allow holons and spinons to hop coherently on to nearest- 
neighbor (nn) sites, which mixes up the two sublattices. 
Different choices lead to different metallic states. The 
state with Bij = B, Dij = D, has spiral magnetic order 
with an incomensurate wavevector. The spiral state is 
energetically favored over other relevant nn states, e.g., 
the flux state, the ferrimagnetic states, etc [23|. Actu- 
ally there is no sign of a spiral (or any other nn state) 
has been found cuprates, or for that matter, in numerical 
treattmcnts of the model. 

Incidentally, the slave-boson RVB theories are also 
based on nn states, with a similar set of order param- 
eters, although interpretations are different [13, [HI- In 
both cases, the nn states are very complex, characterized 
by several order parameters {A,B,D and < 5 >). This 
leads to a complicated phase diagram [l^, [ii| , containing 
several metallic phases, more complicated than that seen 
in cuprates. Furthermore, since slave (or. Schwinger) 
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bosons condense, the physical electron Green's function 
has a pole, with a residue oc | < b > p. Hence the nor- 
mal state is in effect a Fermi liquid at low T, with gapless 
spinon and holon excitations. 

Also, since slave-boson states do not work at half 
filling we should see a transition accompanied by a sta- 
tistical transmutation upon doping. To our knowledge, 
no experimental signature for such a transition has been 
found. Moreover, Haldane and Levine [2^, and inde- 
pendently. Read and Chakraborty [1^ used Berry phase 
arguments to show that, in the presence of short-range 
RVB pairing, spinons are bosons and holons are spinless 
fcrmions, as in our case. Finally, energetics also favor the 
Schwingcr-boson representation, for the simple reason 
that, given the same spectrum, fcrmions will have much 
higher energy than bosons because of the Pauli exclusion 
principle. Therefore, bosonic spinons are better for the 
exchange term. Interestingly, they are also found to be 
better for the hopping term for a; < 1/3 even if J = [sO] ■ 

Instability of spiral states 

Moving holes induce ferromagnetic correlations, which 
compete with antiferromagnctism. In the spiral state the 
system compromises by tilting neighboring spins, so that 
the correlation is partially ferromagnetic. While this al- 
lows the holons to delocalize and reduce kinetic energy, 
there is a substantial cost in exchange energy. Therefore 
such a state may not be stable. Within a Hartree-Fock 
approximation, several authors have found an instabilty 
in the 2 — d Hubbard model toward localization of holes 
in domain walls [3l| for small x. An instabilit y to ward 
phase separation was earlier found by Visscher [32| . 

By comparing the Hartree-Fock energies Hu et al. 
showed that spiral states are indeed unstable relative to 
both phase separation and insulating domain wall states 
in the entire region of the parameter space of interest . 
The actual region of instability is likely to be even larger 
since contributions from RVB singlets to the exchange 
energy are not included in the Hartree-Fock treatment. 
Similar results were found for the t — J model, where 
both antiferromagnetism and singlets are included [27t . 
Since the spiral state has the lowest energy among the 
nn states, it follows that none of these states are stable. 

Longer range Coulomb interaction would destabilize 
the domain walls and spin-charge separated states, and 
at very low doping could stabilze a Wigner crystal of 
holes. However, since t > J the localized states cost 
too much kinetic energy, a metal must emerge at higher 
doping with or without such long-range repulsion. It 
follows that such a metal would have a lower energy 
than spiral state, and can only appear via higher order 
hopping processes by rcnormalized holes, as indicated 
by the one-hole calculations. To derive a rcnormalized 
Hamiltonian and for the results obtained in this paper 
it is not necessary to know the details of the renormal- 
ization procedure. But for the sake of consistency, and 
for a more detailed theory, it is useful to consider the 
possible physical processes involved. 



C. Destruction of Magnetic Order and Renormalization: 

Our work is based on an earlier paper [20| in which a 
self-consistent perturbation expansion in powers of hop- 
ping to one-loop order was used to study the destruction 
of magnetic order upon doping, which is expected to oc- 
cur, but hard to show theoretically. The point is that 
localization of electrons into spins (moments) costs con- 
siderable amount of kinetic energy. As discussed above, 
trying to recover the energy via spin-charge separated 
nn states does not work. Another process which always 
exist is the binding of spin and charge into physical elec- 
trons (or, holes). These try to propagate coherently and 
restore the Fermi liquid, which has low kinetic energy 

In ref. [20[, this process has been treated by RPA, 
in which the electron is treated as a collective excita- 
tion. The zeroth order Green's function Gco is simply 
the convolution of the spinon and holon Green's function. 
Summing the bubbles leads to the full electron Green's 
function 1201 



Gco(kw) 



l-e(k)G,o(k,a;) 



(11) 



where e(k) is the energy of the nonintcracting electron. 
This can be viewed as a generalization of the one-hole 
method where a similar expansion is carried out for one 
hole. In the many hole case, the effect of the electron 
on the spinon and holon self-energies is also calculated, 
using Gc- Unlike moving holons which only change the 
direction of the local moments, a moving electron has a 
much more violent effect, as it also tends to delocalize, 
and hence destroy the moments. It was found that x > 
Xc, where Xc is rather small, long-range magnetic order 
is destroyed at T = 0, i.e., there is no Bose condensation 
of spinons [13] . 

It was also found that Gc does not have a pole at low 
frequencies, that is, the system is not a Fermi liquid for 
small X. The reason is that in order to move an elec- 
tron must avoid other electrons of opposite spin (30j . For 
small X and at low dimensionality, such pathways are 
rare. Consequently, the electrons (or, physical holes) do 
not become coherent. However, spinons become gapped 
since long-range magnetic order is destroyed. The gap 
shows up in Gco and hence in Gc [33j . In contrast, in the 
one-hole case, the Green's function has a pole, which is 
not surprising since a single hole can not destroy long- 
range AF order. In either case, the holes are strongly 
rcnormalized by the short-range incoherent hops. 

Interestingly, it was also found that the presence of a 
RVB condensate leads to a second order hopping process 
in which a pair of holons hop together, accompanied by 
a singlet backflow, which automatically leads to super- 
conductivity. This is the pairing mechanism considered 
below. The theory was further developed in ref.plj 
where an effective Hamiltonian for the rcnormalized 
holons including the pairing term was derived. However, 
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it was not entirely correct - since the normal state was 
assumed to be a spiral metal, with nonzero B and D, 
but without long-range magnetic order. Such a state is 
not stable. Below we give a derivation which corrects 
this flaw. 

IV. Renormalized Hamiltonian for Small Doping 

The spinon gap allows us to obtain an effective Hamil- 
tonian involving renormalized holes and singlets, by in- 
tegrating out the spinous and setting w = as was done 
in refs. [20l. [2H. Here we give a simpler derivation using 
the usual continuity arguments Since the nn MF 

state is not stable we have B = Q and D = 0. A moving 
hole affects antiferromagnetic configurations and singlets 
differently. In the former string of ferromagnetic 

bonds are created whose energy increases with the length 
of the path, and which can not be easily repaired quan- 
tum mechanically. The hole is essentially forced to re- 
turn, creating a renormalized hole localized within a bag 
of spin excitations. This process have been considered 
by many authors in the context of the onc-holc problem 
[24 ,25i , ,26j . ft is also the reason for localization of holes in 
domain walls in the Hartree-Fock treatment of the Hub- 
bard model [1^, [3l| . The hole moves rapidly within the 
bag, gaining kinetic energy of order t. This would con- 
tribute an incoherent (essentally, k-independent) back- 
ground to the electron spectral function as well as to the 
optical conductivity. We assume that disordering of the 
spins prevents long-range AF order, and leads to short- 
ranged AF correlations. The renormalized hole is carries 
the spin excitataions, and is thus larger and more mas- 
sive. 

After this intital renormalization, by continuity the 
Hamiltonian would be similar to the original one, except 
for the absence of magnetic order, as long as no addi- 
tional symmetries are broken, ft will be characterized by 
a hopping amplitude <e// < ^, and an exchange coupling 
Jeff- For small x, we expect Jeff ~ J- Also, one-hole 
calculations [H, , suggest that tef / scales with J for 
t > J. Such estimates can be used since the renormaliza- 
tion process involoves only nearby hops, and therefore at 
low-hole densities it is essentially a one-hole problem. An 
important aspect of the one-hole calculations, is that the 
local constraints are imposed exactly. The constraints 
obeyed by the renormalized particles are thus expected 
to be weaker, although the renormalized Hamiltonian re- 
tains gauge invariance. 

The singlets are affected less drastically by hole mo- 
tion since the offending configuarations can be repaired 
more easily, as discussed below. We therefore assume 
that the problem is already renormalized by the fast pro- 
cesses described above. When a renormalized hole hops 
from sublattice a to 6 it breaks a singlet and creates two 
spinous, costing, say, an energy > 2As. The system 
relaxes by a second hop, after which the singlet is re- 
constructed. In the low frequency limit {uj << f2) these 



processes can be described by 

^2 

H,nt = -^<V J2 cLc..cL,c,.,P >, (12) 

where Cia = bl^hi] but now all the operators correspond 
to renormalized particles. Here < V.-.V > means that we 
keep only those terms that involve renormalized singlets 
and holes. Although nominally similar to superexchange, 
the intermediate state does not involve double occupancy 
which has already been projected out. As shown below 
the new interaction terms can be rewritten in terms of 
singlet operators Ajj . Now, Aij and Bij decay exponen- 
tially as e"'''^/^^ [13], where ^ is the magnetic correlation 
length. Therefore, to obtain a minimal Hamiltonian, it 
is sufficient to retain only the short-range hopping terms 
since longer range terms are exponentially suppressed. 
The longer-range hopping terms also preserve the under- 
lying symmetries, and neglecting them should not change 
the physics qualitatively. 

Once the singlet is broken, there arc three ways to re- 
move the excess energy. First, the hole can hop back, and 
the singlet reconstructed. This is confining and its effect 
is to renormalize the chemical potentials. The remaining 
two processes can lead to hole propagation. 
A. One-hole process: The hole hops to another site on 
sublattice a, and the singlet is reconstructed on a differ- 
ent link [Fig. 1] . Notice that this process involves succes- 
sive hopping by two electrons of opposite spin, but only 
one holon. Projecting onto the spinon singlet subspace 
[P{bl^bj^^cr'P = crAy, etc], we obtain 

ijl 

where, A^ creates a singlet and h'^ creates a renormalized 
holon, and ts = is the effective intra-sublattice hop- 
ping parameter. Now, the term (1 — h^hj) arises from the 
intermediate excited state. It is gauge invariant, and for 
small X, can be approximated by its average value 1 — x, 
leading to 

^^-^{l-x)Y,AlA.,,hlhi, (14) 

ijl 

This term describes hopping by a holon within the same 
sublattice, accompanied by a backflow of spinon singlets. 

The one-hole term above differs from the usual next- 
nearest- neighbor hopping {t'). In the early RVB-gauge 
theories based on similar RG considerations it was as- 
sumed that the result of second order hopping is to gen- 
erate an effective t' term of the form 

t'^cla^^la +h.C., 

a 

where i and I are next nearest neighbors. To generate 
such a term from Eq. (12) would require consecutive 
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; i \ o o 

FIG. 1: One hole process. A hole hops from I to j, breaking 
the singlet (ij) denoted by the solid line. Here i and I are on 
sublattice a and j is on sublattice b. The hole then hops to i, 
and the singlet is reconstituted at (jl) 

hops by two electrons of the same spin. This is not eas- 
ily done when the initial and final states belong to the 
singlet subspace, as in our case. For after the first hop 
the remaining electron will have opposite spin. Then to 
generate a t' one has to flip both spins via additional 
higher order intermediate processes, which arc not very 
likely. (It can be done via J' interaction. But for this 
model J' is usually much smaller than J) . Of course the 
actual model that applies to real cuprates will have a t' 
term to start with which, as will see later, will play a role 
since it also preserves the two-sublattice character, but 
it would not replace the dominant one-hole term derived 
above. 

B. Two-hole process: The system also relaxes if a sec- 
ond hole hops from sublattice b to a, and the singlet is 
reconstituted on a different link [Fig 2]. This process 
yields a term 

~tsJ2AliM,hlhlhih^, (15) 

ij:lrn 

which describes hopping by a singlet, accompanied by 
the backflow of a holon pair. This is the small x form of 
the interaction derived earlier plj . but here the normal 
state is different because of the one-hole term above. 

i O' o t o 

FIG. 2: Two hole process. A hole hops from I to j, breaking 
the singlet (ij) , a second hole hops from m to i, and the singlet 
is reconstructed at (ml). 

The full Hamiltonian also includes the rcnormalizcd 
exchange term, characterized by Jeff- For the pur- 
poses of this paper it is sufficient to treat Jeff, teff 
and 17 as phenomcnological parameters. However, as 
discussed earlier, one-hole calculations [1^1 suggest that 
they all scale with J. Hence the renormalized Hamilto- 
nian describes a strong-coupling problem. For cuprates, 
J ^ 0.\eV, and is the largest energy scale in the renor- 
malized problem. If numerically t(.ff > then a fre- 
quency cut-off can be used. By continuity qualitative 
physics will not be changed. As it turns out (see below). 



actually the effective holon hopping amplitude teff A is 
~ JA T* . Since A is only a small fraction of unity and 
decreases with x, a cut off may not be necessary, at least 
for some values of x. 

Because of the spin gap our renormalized Hamiltonian 
describes a short-range RVB model for small x p^ . 
However, the SR Hamiltonians used by earlier authors 
[l9| were different (for example, the second term which 
is responsible superconductivity was absent) since they 
were not required to agree with Arovas-Auerbach state 
at half filling. Also, holons were thought to be bosons. 
It was later shown that for short-range RVB models 
they are actually fermions, which agrees with our choice 
[28l . . Our renormalized Hamiltonian can not possess 
a conventional Fermi liquid state since the latter has 
gapless spin (and charge) excitations. 

V. Analysis of the Renormalized Hamiltonian 

A. Non-ordered Phase: Strange Metal or Insulator?: 
There is always a state in which no symmetry is bro- 
ken (even at the MF level). Then singlets are not con- 
densed, i.e., < Aij >= 0. Hence, 

< >= 

for i 7^ j, so that there is no coherent propagation 
of holons. Holons are localized due to strong gauge 
fluctuations. This state has no long-range order, and 
would always occur at sufficiently high T (above T* 
in our case). We identify it with the so-called strange 
metal phase of the cuprates. Since electrons are also 
incoherent, there are no coherent charge carriers at all in 
this state. Therefore, the strange metal is not a metal; it 
is a new type of quantum insulator. In effect, this state 
corresponds to a Gutzwillcr projected Fermi sea; but 
one with strong singlet fluctuations, which lead to the 
gapping of the spinous (and, hence, the electrons), and 
non-Fermi liquid behavior. In contrast, in other RVB 
theories the state above T* is metallic, with gapless spin 
and charge excitation [l^, [ll|. In these theories the 
nonordered state appears at a much higher temperature, 
and is ignored. 

B. Pseudogap Phase 

To proceed further, we first do a MF decomposition 
of Hint to obtain two separate Hamiltonians for spinous 
and holons, which breaks gauge symmetry. Then, from 
the one-hole term, we get 

ijl 

(16) 

In principle, we can find a state for which < Aij >= 
(no singlet condensation), but < aLa^ >^ 0. Such a 
state will not be favored energetically, since the energy 
gained 0{tsx) per bond is low compared with J^ffA^ 
gained from singlet condensation. 
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This leaves the RVB state which becomes stable below 
T* with < Aij >^ 0. From Eq. (2) we obtain 

Using this in (15) we obtain an effective holon hopping 
Hamiltonian 

Hha = - ^ th,ijh\hj, (17) 

which describes coherent holon propagation within the 
same sublattice. Here th^ij — ±th for i,j next-nearest 
neighbors along the (1,±1) diagonal direction, t^^ij = 
t/i/2 for next-next-nearest neighbors (along the x and y 
directions), and = tsA^{l — x). The holon energy is 
then given by 

eh.i(k) = -2th + 2th{sink^ + smkyf. (18) 

The holon band (hence, metallic conduction) appears as 
soon as A ^ without additional symmetry breaking. As 
expected, the spectrum has the two-sublattice character 
since e/ii(k) = ehiiQ - k). 

We can do a similar MF decomposition of the two-hole 
term (Eq. 14), and replace the Ay's be their expectation 
values, which yields a holon- holon interaction term of the 
form 

Hh,^nt = -tsA^ Py-F™/ + '^■^1 (W) 

ij:lrn 

where fJ^ = creates a holon pair on the link ij, 

and the sum is over plaquettes. The order of the indices 
follows from the symmetry of Aij and is very important. 

Once the holons are propagating, the pair-hopping 
term will also contribute to the normal state holon spec- 
trum. Doing a MF decomposition of Hh,int we obtain a 
term 

H2 = t,A^ [Duh]h^ + D^jh\hi - DuD^.j + h.c], 

ij,ml 

where 

Du =< h\hi >= ^ E e^"^-^"^'"""' V(e/> (k) - m), (20) 

k 

is the average holon hopping amplitude, and / as the 
Fermi-Dirac function. Using the symmetry of the holon 
spectrum, we find that Du = 0, when I and i are nearest- 
neighbors. And, Du = ±Di for I and i along the (1, ±1) 
diagonal directions, respectively, with 

Di = --^ X!^"^^^ sinfc.y/(eft(k) - /i). (21) 

k 

Therefore the extra diagonal hopping terms have the 
same symmetry as the original ones. Furthermore, since 
the minima of the holon spectrum is at ±^{11, — tt), Di is 



positive. Using the fact that hopping along (1,1) is ac- 
companied by a backflow along (1, —1) the extra diagonal 
terms are found have the same signs as before. Including 
these the full holon spectrum becomes 

e/i = — 2t/i + 2t/i(sin + sin fcy)^ -I- 4Dits sin k^ sin ky. 

(22) 

The parameter Di has to be calculated self-consistently. 
However, for small x, the main contribution comes from 
near the minima of the spectrum ( k = ±(7r/2, — 7r/2)), 
where we can set sin kx sin ky = —1, which yields Di w x. 

So far we have considered only the holon part. The 
motion of holons (or holon pairs) is accompanied by a 
backflow of singlets. The mean-field singlet Hamiltonian 
is then given by Hh = Hj + Hbi + Hb2, where Hj is 
the usual exchange term (see Eq. 1), with an exchange 
constant Jeff] 

Hbi = - a:) E AtjA,, < hlhi >, (23) 

iji 

is the singlet backflow term associated with the one-hole 
process, and 

Hb2 = -ts Y ^L/Ay < h\h]hiK, >, (24) 

is that associated with the two-hole process. The last 
two terms only appear below T*. They will contribute 
to the singlet condensation energy in the pseudogap 
phase. They are suggestive of the emergence of coherent 
singlet excitations, which are spin-0, chargeless vector 
bosons. These will contribute to low-T specific heat, and 
thermal conductivity, but not to electronic transport. 
We will study them in a future paper. 

Symmetry Properties of the Pseudogap Phase 

The pseudogap state preserves the two-sublattice (red- 
blue) character of the RVB state. To see this consider 
the correlation functions for the renormalized particles. 
Using e/j(k) = e/j(Q - k), we find Aj =< h\hi >= 
cos(Q.(rj — Tj)), which is thus nonzero only on the 
same sublattice (even). The spinon correlation functions, 
and hence the gauge- invariant magnetic correlation func- 
tions have the same symmetry as in the Mott phase. 
The electron hopping amplitude Pij^a =< c\^Cja >= 
—BijDij. Since Bij and are both even, is nonzero 
only on the same sublattice. The Fourier transform of Py 
is the momentum distribution function, which thus sat- 
isfies nc(k) = ric(Q ~ k). In a Fermi liquid the electron 
hopping amplitude shows a power law decay. In our case, 
the holon hopping amplitude Dij shows a metallic, i.e, 
power-law behavior since holons are gapless. However, 
since spinous are gapped Bij, the full electron hopping 
amplitude Pij , decays exponentially, refiecting non- Fermi 
liquid behavior. 

The metallic character of the holons can also be 
seen in the gauge-invariant charge structure factor. 
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Let Pi = hjhi— < hlhi > measure the excess hole 
density. Then the charge structure factor is given 
by: Sch.ij =< PiPj >= -lAjP, for i ^ j, and 
Sch,ii = 2^(1 — x). This is nonzero on the same 
sublattice, and exhibits the long-range oscillatory struc- 
ture of a metal. In the momentum space it satisfies 
5c/i(k) = 5c/i(Q - k) (Fig. 3). In contrast, 5c/i(k) 
of an ordinary metal increases from zero at k = 
and becomes a constant for q > 2kp. In our case, an 
image of the behavior near k = appears near k = Q. 
Experiments probe the bare correlation functions. These 




FIG. 3: Charge structure factor, showing thetwo-sublattice 
symmetry. The calculation is approximate. Note the symme- 
try between k = and k = (tt, tt), which is distributed among 
the four corners of the Brillouin zone. The valley region is el- 
evated relative to k = 0. In a normal metal, the structure 
near (7r,7r) is absent. 

are dominated by short-range incoherent processes that 
do not preserve the two-sublattice property, which is 
thcrfore not easy to see. The best candidate is iSc/i(k) 
since holon motion is coherent. The experimental SchO^) 
would no longer vanish at Q, but there will still be a 
dip. 

VI. Superconducting State 

It is clear from the structure of the two-hole term (Eq. 
19) that the system can lower its kinetic energy if two 
holons hop as a pair. Pairs will condense, leading to 
Fij =< Fij 0. Let us define the pairing order pa- 
rameter for physical electrons as a singlet 

(cjic^t - CjTQi)/2 > . (25) 

Then, Cij = —AijF*j ^ since the spinous are already 
condensed, giving rise to superconductivity below Tc < 
T*. The order parameters Aij and Fij arc not gauge 
invariant, but Cy- is. The symmetries of Aij and the 
holon spectrum are of course already known. 

We solve the superconductivity problem by a mean- 
field approximation. Since holons arc fcrmions the holon 



pairing order parameter satisfies Fjm = —Fmj- We de- 
note the vector Fjm by two components: Fjm ~ iFx if Tn 
is to the right of j, and Fjm = iFy if m is above j in the y 
direction. The prefactor i is chosen so that order param- 
eter for the electron (Cooper) pair is real. For a uniform 
system, we can take \Fjm\ — Fq, but the phases along x 
and y need not be the same. Without loss of generality 
we can choose, F^ = Fq and Fy = aFo, with a = e*^. 
The choice of a = ±1 leads to Cx = ^Cy, corresponding 
to s-wave (d-wave) symmetry for the electron pair wave 
function. But other values of a can, in principle, lead to 
s + id OT s + d symmetries. To find the correct symmetry 
we compare the free energies for different choices. The 
MF Hamiltonian is given by 

Hmf = J2 iC^K^k + ^Ah{k){hlhl^ + h.c), (26) 

k 

where ^(k) = (k) — /i, p is the holon chemical potential, 
and A/i(k) — 2tsF(k) is the holon gap function, with 

F(k) = 2FQ{si\ikx +as\nky). (27) 

The holon spectrum e;i(k) is assumed to contain the 
Hartree contribution from the interaction term. The 
Hamiltonian is diagonalized by the Bogulyubov trans- 
formation. The self-consistent equations for the order 
parameters are then given by 

. , , . , . , ,tanh(/3i;(k)/2) 
F^ = — 2_^smfc^(smfc3: + asmky) -^^^ , 

(28) 

where rj = [x, y), and 

i?(k) = [C2(k) + |A,(k)p]V2. 

Since F^ is real wc find that solutions exist only for 
smO = 0, i.e, for a ~ ±1. Combining the equations 
for Fx and Fy, we obtain 

1 Iv-T,.., , ■ , ,2tanh(/3£:(k)/2) 

(29) 

VF(k) is a suitably chosen cut-off function. 
To find the symmetry let us consider the T = case, for 
which tanh/3i?(k)/2) = 1. We have solved this equation 
numerically with and without a cut-off function. Quite 
generally, we find that a = — 1, (i.e., d-wave) is favored, 
as it leads to the largest Fq, and hence the largest con- 
densation energy. The origin of this result can be un- 
derstood from the following simple considerations. The 
dominant contribution to the sum comes from the region 
where |^(k)| = |eh(k) — p.h\ is small, and the symmetry 
factor I sin kx + a sin ky \ is large. As shown in Fig. 4, the 
holon Fermi surface is in the second and fourth quadrant, 
exactly where sin kx+a sin ky has maxima for a = —I {d- 
wave) and vanishes for a — I (s-wave)). Hence, d-wave 
always wins. Thus the symmetry of the superconducting 



11 



3 




-3- 

-3-2-10 1 2 3 
kx 




FIG. 4: Origin of d-wave symmetry. Upper panel: Holon 
Fermi surface. Holons live within the pockets centered at 
^(tt, — tt) and i(— 7r,7r). The shape of the FS is somewhat 
different when compared with ref. because the holon spec- 
trum has corrections from the two-hole term. The FS is not 
directly observable since under a gauge transformation it is 
moved and deformed. Lower panel: The symmetry factor 
I sin kx + a sin ky\. It has broad maxima in the hole-rich re- 
gion for a = —1, resulting in ma^cimum condensation energy. 
In contrast, the symmetry factor vanishes in this region for 
a = 1 (s-wave). Under a gauge transformation both the Fermi 
surface and the symmetry factor move together to preserve 
these gauge-invariant results. 



order parameter is determined by the symmetry of the 
underlying RBV state. 

Note that since F(k) = i^(Q — k), the two-sublattice 
property is also preserved in the superconducting state. 
The holon pairing function is odd since for any two i,j 
satisfies Fij — —Fij cos Q.(ri — rj). Since Aij is also odd, 
the electron pairing function Cij = —AijF*j is also odd. 
Similarly, the symmetries of nc(k) and spin-spin correla- 
tion function remain unchanged in the superconducting 
state. The charge structure factor, however, picks up an 
additional contribution: Sch,ij = l^y P ~ l^yP: s-^d is 
no longer restricted to the same sublattice; but like the 
spin-spin correlation function, it oscillates in sign. 



Note, the rcnormalizcd Hamiltonian describes a 
strong-coupling problem since the onc-holc and two-hole 
terms are essentially governed by the same energy 
scale tg. The holon wave functions depend on x; but 
not on ts- Hence the BCS like MF approximation is 
expected to work only for T « Tc- In particular, 
treatment of the transition region, and the calculation of 
Tc{x) will require strong-coupling techniques. Similarly, 
MF approximations would not work near T*, where 
fluctuations will be strong. The phase diagram shown 
in Fig. (5) is thus only schematic. The main point 
is, unlike earlier RVB theories, the number and type 
phases for small x are similar to those in cuprates. 
Moreover, as discussed below, the qualitative behavior 
of various phases can be easily understood. Given 
the numerous unusual properties of the cuprates, and 
the highly constrained nature of the rcnormalizcd 
Hamiltonian, the predictions can be put to severe ex- 
perimental tests. So far we have found no contradictions. 
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FIG. 5: Schematic phase diagram. The state above T* is 
a novel insulator and is unordered. The pseudogap metallic 
state is a spinless Fermi liquid. Except for the strange insu- 
lator, the remaining phases are characterized by RVB order 
parameter A^Q 

VII. Comparison with Experiments: 

Clearly, the present theory is consistent with the 
important features of the experimental phase disgram. 
By construction, it is consistent with the physics at half 
filling. The normal state has spin-charge separation and 
is not a Fermi liquid. There is a transition from a strange 
phase to a pseudogap metal below T* . As discussed 
later in section VIII, these phases are two-dimensional, 
whereas the superconducting state exhibits 3d behavior. 
However, the phases themselves are predicted to have 
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novel chracterisitcs that defy conventional wisdom. The 
key ones are: (1) there is a true spin gap O which 
exists in all three phases, and which is distinct from 
the pseudogap that appears below T*, (2) holons are 
confined above T*, so that the strange phase is an 
insulator rather than a metal; (3) below T* holons 
form a spinless Fermi liquid of concentration without 
an observable (small) Fermi surface. The pseudogap 
metal and the supercondor arc characterized by the 
RVB order parameter and the two-sublattice property, 
whereas the strange phase does not have any order. 
These results are robust as they follow from symmetry 
and particle statistics. A large number of important 
conclusions can be drawn from the very structure of the 
effective Hamiltonian in each phase, which do not require 
detailed calculations, and which can be tested against 
experiments. Wc now examine recent experimetal data 
to show that there is very strong evidence in support of 
these predictions. We separate the experimental findings 
into three groups: the spin sector, the charge sector and 
the electron sector. 

A. The Spin Sector: 

The spin gap implies the existence of an additional 
temperature scale ~ n/k. In the pseudogap region 
the AF correlation length ^spin is expected to be short 
- no more than a few lattice spacings - due to strong 
renormalization processes implicit in our theory. Then, 
n ~ Jeff{x). On the other hand, T* - Jeff A < T", 
since A is only a fraction of unity (from the MF the- 
ory). At half filling, RVB ordering arises from long-range 
singlet-singlet interactions mediated by gapless spinous. 
In the doped region this interaction is weaker and short- 
ranged because of the gap. However, now the renor- 
malized hopping terms favor singlet condensation since 
holons can propagate coherently only if singlets are cor- 
related. The qualitative behavior is easy to understand. 

Above we expect free spins to dominate. Then 
the uniform paramagnetic susceptibility Xpara should be 
Curie like, i.e., decrease with increasing T. Below 
singlets will form, and Xpara will start decreasing with 
decreasing T. Note that there is no transition at T^] 
it is just a broad crossover scale. Hence, we can asso- 
ciate it with the maximum of Xpara- With decreasing T, 
there will be more singlets and they will become increas- 
ingly more correlated, and eventually condense at T* . 
Below T*, Xpara will dccrcasc much more rapidly, and 
vanish as T — > (because of the gap), even in the ab- 
sence of superconductivity. The magnetic behavior above 
T* is determined by fluctuating singlets correlated over 
a distance ^singlets (T) , together with weakly correlated 
free spins. The behavior near T'^ and above, where sin- 
glets are weakly correlated, can be crudely understood in 
terms of a one singlet (two-site) problem, for which Xpara 
is given by 

Xpara2 - 3 ^ ^^j^^^ ■ 



This has gap fl = Jeff-, and it vanishes exponentially as 
e~0J^ff as r ^ 0, shows a Curie like 1/T behavior at 
large T with a maximum at T'^ « 0.62 Jef f /k. Fig. 6 
shows Xpara2 SlS a fuuctiou of T. 




FIG. 6: Paramagnetic susceptibility for a two-site problem. 
Note the maximum at T*^, the spin-gap scale. The line on the 
right below T* is schematic, showing how singlet condensation 
at T* for a many-site problem would further depress Xpara- 

The pseudogap has been observed in cuprates by using 
nuclear magnetic resonanancc (nmr) techniques [3, Issj . 
There is a rapid decline in Xpara, as measured by the 
Knight shift, below T*, which is far above Tc- A similar 
downturn is seen in the spin-lattice relaxation rate, 1 /Ti . 
More interestingly, the higher crossover temperature 
was also found in many of these materials (for a review, 
see 0). That T° has not been seen in some materials is 
more a matter of not going to high enough temperatures. 
For example, was not seen in initial measurements on 
YBa2Cu408 upto 400 K. Curro et al.^ extended the 
measurement to 715K and found a broad maximum at 
T° ^ 500i^, with Xpara decreasing slowly on both sides 
of T". In this material T* is about 240X, and Tc - 81A'. 
It is interesting to note that To and T* vary roughly the 
same way, decreasing with increasing x, consistent with 
the theoretical expectation that they arc proportional to 
the same energy scale Jeff{x). 

Our theory also predicts that Xpara would vanish 
in the normal state as T ^ 0. This is harder to 
confirm since supercondcutivity intervenes. However, as 
discussed in ref. 0, in highly unde rdop ed systems such 
as YBC06.7 [H and YBa2Cu408 [13, there is almost 
no sign of a superconducting transition in the nmr data 
at Tc, indicating that there is no additional spin pairing 
at Tc, although Xpara continues to decrease. Similar 
results have also been seen in moderately underdoped 
Bil222 [3^ by Walstedt et al. who argue that the lack 
of any effect at Tc indicates that the gap is unrelated to 
superconductivity, and represents spin-charge separation. 

B. Electron Sector: Tunnehng DOS 
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Further support for the theory comes from tunneling 
experiments, which probe the electron spectra (charge 
plus spin). Unlike magnetic properties, tunneling DOS 
j4ll . |43 | shows a sizable effect near Tc in underdoped 
Bi2212. For T « there is a gap in the DOS, with 
well-developed conductance peaks at ±Asc, except that 
2Asc/kTc is much larger than the BCS value. As T ap- 
proaches Tc from below, conductance peaks decrease in 
height and disappear at Tc, as does the zero-bias peak 
originating from the Josephson current. This is what is 
expected. However, while the gap decreases as is ap- 
proached, it does not close; instead a gap like structure 
with a depressed DOS continues to exist above Tc- Al- 
though, sometimes it is referred to as the pseudogap in 
the tunneling literature, the gap actually exists far above 
T* (as measured in the nmr experiments), that is, it also 
exists in the strange phase up to high temperatures, i.e., 
it is the spin gap. Such a behavior can not be understood 
in terms of nn (e.g., slave-boson states, or a Fermi liquid. 

The tunneling data are clearly consistent with our the- 
ory in which the spin gap exists in all three phases. The 
gap and the singlet condensation qualitatively account 
for the nmr data in the normal state. Superconductivity 
arises from a pairing of spinless holons. Since there is 
spin-charge separation Xpara 

and 1/Ti would largely be 
unaffected below Tc, as seen in the experiments. How- 
ever, tunneling experiments measure the DOS of the elec- 
tron, containing both spin and charge. At the MF level 
this is a convolution of spinon and holon DOS. Above 
Tc holons are gapless, so the gap in the electron Green's 
function reflects the spin gap. This is why tunneling gap 
exists far above T*. Below Tc, a gap also opens up in 
the charge spectrum, as holons pair; this explains the in- 
crease in the total gap size with decreasing T, and the 
appearance of the conductance peaks. Also, experimen- 
tally the T = gap increases with decreasing a;, mirroring 
the behavior of T° associated with the much larger spin 
gap, as expected from the theory. 

Experimental observation of a spin gap far above T* 
rules out a mechanism based on preexisting electron 
pairs, carrying both spin and charge. To account for 
the gap, such pairs should be formed near T", and pre- 
sumably condense below Tc- It is then hard to explain 
T*. There are other difficulties. The condensation at Tc 
must increase the binding energy of the Cooper pairs, 
enough to cause the large depletion observed in the DOS 
below Tc, and the observed increase in the size of the 
gap. However, since electrons forming a Cooper pair also 
carry spin, such a large increase will be accompanied by 
a large change in Xpara and 1/Ti below Tc, a behavior 
not seen in the experiments [35| . 

C. The Charge Sector 

Strange Insulator: A key prediction of our theory is 
that holons are confined (nonpropagating) above T* since 
A = 0. Furthermore, since electrons are incoherent and 
gapped the 'strange metal' is not a metal, but a new type 



of insulator. A strong evidence for this comes from the 
measured in-plane dc resistivity, p. In a metal p is sup- 
posed to saturate at the Mott value. Actually, while the 
linear T dependence of the p{T) has received most of the 
attention, it has been pointed out by many authors that 
p shows no sign of saturation, and in the underdoped 
regime, far exceeds the Mott value psi. [4^. l45j . The fact 
that p{T) increases with T does not make the system a 
metal since p also increases in a disorder-induced (An- 
derson) localized insulator as long as inelastic mean-free 
path iineiiT) is less than the localization length, ^loc- 
A quantitative theory is beyond the scope of this pa- 
per; but the qualitative behavior can be understood, as 
follows. In the case of Anderson localization ^loc is inde- 
pendent of T; then p{T) would show an upturn at low 
T, as iinei{T) exceeds ^loc- In our case, the appropri- 
ate localization length is the distance over which holons 
can move coherently, which is the singlet-singlet corre- 
lation length £,singiet (T) (not the shorter AF correlation 
^spiri)- Since £,singiet(T) incrcases with decreasing T, and 
becomes infinite at T* (at the MF level) where £inei is 
finite. Hence there can not be a low-T upturn, which is 
consistent with experiments. On the other hand, since 
there is no coherent holon motion in the perpendicular 
direction the c-axis resisitivity should show an upturn, as 
seen. 

A clear indication is also obtained from the frequency- 
dependent conductivity (7{uj) which, for a metal, should 
exhibit a Drude peak. But in our case there would be no 
such peak above T*. Recent experiments in LSCO (isj 
and B21212 [i^ show that such a peak does not exist 
for small x. Instead, one has nearly a flat spectrum over 
a range of frequencies of order 1 cV. The lack of a Drude 
peak indicates the absence of coherent charge carriers. 
The flat spectrum would come from incoherent hopping. 
At higher doping, a broad peak develops, which is not 
surprising since in this case the physical electron comes 
into play as one approaches the Fermi liquid. 

Pseudogap State: Emergence of Metallic Conduction 

The theory predicts that below T* holons become co- 
herent and form a spinless Fermi liquid of concentration 
X. Therefore p(T) should drop rapidly with decreasing 
T, and become metallic. Indeed, ex per imentally p{T) 
is found to drop faster than linearly [43l . IZtI . [49| below 
a temperature which approximately agrees with T* ob- 
tained from nmr experiments 

U,m- At lower T, p{T) 
becomes metallic and has a residual impurity contribu- 
tion dl, m, [ill, as expected of a metal. At low-T, 
p{T) oc T^ due to fermion-fermion scattering [H, Isot . 
also as expected. For LSCO one also finds that p does 
show an upturn at very low T, which is believed to be 
due to disorder- induced localization [isj . 

Recent experiments [10, [Ml in underdoped LSCO and 
YBCO show that the Hall coefficient Rh is independent 
of T and proportional to 1/x in the pseudogap regime, 
as expected of a metal. This should be compared with 
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the strong T-dcpendence in the 'strange' insulator phase. 
With increasing temperature, p{T) is seen to deviate 
from the to eventually a linear-T behavior above T* . 

Strong evidence also comes from optical conductivity 
cr{uj), which shows a Drude peak at low J\ with an in- 
tegrated area (spectral weight) oc x It is char- 
acterized by a small plasma frequency, consistent with 
a small bandwidth, or heavy holons, as in our theory. 
With increasing T, a{uj) broadens, but at a rate too 
large to be attributed to thermal effects in a Fermi liq- 
uid, and merges into the incoherent background above 

These transport properties, taken together, strongly 
support the view that the charge carriers are holes, 
rather than electrons, and they form a Fermi liquid 
below T*, and when combined with the nmr data, they 
suggest that these holes are spinless, supporting our 
view. They also suggest that at small x there are no 
coherent charge carriers above T*. 

Lack of a Fermi Surface: 

If the pseudogap state is a Fermi liquid of concentra- 
tion X as the transport experiments suggest, where is the 
corresponding, presumably, small Fermi surface? This is 
a real puzzle. Attempts to understand this (50l | in terms 
of the 'Fermi arc' found in the photoemission experiments 
[ssl . [5^ does not make sense for many reasons. First, 
there are no quasiparticles associated with the Fermi arc. 
The electron Green's function is completely incoherent in 
the normal state; there are no sharp peaks in the pho- 
toemission spectrum. Therefore, the Fermi liquid is not 
formed by electrons. Furthermore, as discussed earlier, 
neither the mmr data nor the gap in the tunneling spec- 
tra is consistent with an electron Fermi liquid. 

In contrast, our theory predicts that the holon Fermi 
surface is unobservable, even in principle. This is because 
the holon Green's function, and thus the holon spectrum, 
are not gauge invariant. Hence the holon Fermi surface 
can be deformed and moved around by a gauge trans- 
formation. The holon Fermi surface does not exist when 
averaged over all the gauge-equivalent copies. This is 
why experiments do not see it. 

Another prediction is that since Xpara decreases 
rapidly below T*, and eventually vanishes as T — > 0, the 
total uniform magnetic susceptibility of pseudogap metal 
would become more diamagnetic with decreasing T due 
to the orbital motion of holons. The effect should be 
small near T* , but should become more prominent in the 
metallic region (at lower T). Unfortunately, supercon- 
ductivity intervenes. Note that the diamagnetic response 
discussed here should not be confused with strongly T- 
dependent response found by Wang {et al. [11] consid- 
erably above Tc- The latter is presumably due to pair- 
ing above T^- Such strong-coupling effects are not ruled 
out by our theory, but its treatment would require more 
sophisticated techniques. We point out that in the an- 
laysis of Wang et al. the normal-state diamagnetic con- 
tribution appear to have been subtracted out along with 



other weakly temperature-dependent terms. Therefore 
more experimental work will be needed to see the effect 
in the normal state. 

In our theory superconductivity appears via a pairing 
of holons. Since both terms in the effective Hamiltonian 
arises from hopping, condensation energy is oc ~tj,/th, 
which causes a reduction in the kinetic energy, as 
observed [56l . [57| . Our mean- field treatment implies 
that spin-charge separation continues to exist in the 
superconductor. For small x, there is some evidence for 
this since Xpara is essentially unaffected at T^, but more 
work is needed to clarify this issue. 

VIII. Origin of Tv^ro Dimensionality 

In this section wc consider the origin of two- 
dimensionality of the normal state(s). which is central to 
the occurrence of high-Tc superconductivity in cuprates. 
These are highly anisotropic 3d materials. Ordinarily the 
metallic state would show 3d behavior below some Ti 
since confinement within a plane would cost too much 
kinetic energy. For cuprates, this would imply T±_ < Tc, 
which is small at low doping, and actually vanishes at 
some critical x. In other words, the normal state appears 
to remain two-dimensional as T ^ (in the absence of 
superconductivity). On the other hand, the other phases 
in cuprates: the AF state at half filling, the supercon- 
ducting state, and the Fermi liquid state at large x are 
all three dimensional, as one would expect. 

In theoretical studies, the 2-d nature of the normal 
state is usually assumed^ not established. One simply 
analyzes the 2d Hamiltonian. The question arises: does 
the corresponding normal state remain two-dimensional 
when hopping is turned on in the perpendicular direc- 
tion? If it does not, the theory should be discarded. 
This requirement, seldom tested, provides a powerful 
constraint on any theory. 

Out of plane hopping can be modeled by adding 

l.Z 

to the original model, where t±_ << t; and for simplicity 
we consider a tetragonal lattice. Suppose we put [/ = 0, 
and treat the anistropic hopping model as an effective 
model for the Fermi liquid state. Obviously, this would 
show 3d behavior. The same is true if the particles are 
bosons, as in a superconductor or a quantum magnet. 
Similary, the nn (e.g., the spiral or the slave-boson) states 
should be three-dimensional since they arc modeled by 
similar effective MF Hamiltonians. in which fcrmions and 
bosons hop 'freely' in all three directions. The issue of 
confinement has been studied for coupled Id Luttinger 
chains which have gapless excitations ^58J . 

In our theory the gapping of the spin-excitations pro- 
vides a protective mechanism for 2d confinement. This 
is far from obvious since, as we have seen, holons can 
delocalize in the plane even though spinous arc gapped. 
So why not in the z direction? First, note that a nonzero 



15 



t± leads to an exchange interaction J± = At\/U = 
J{t±/t)'^ << J. Since a spin can belong to only one 
singlet at a given time it follows that the spin singlets 
are formed only on the x-y plane in the ground state. In 
order to form them along the z direction one has to break 
singlets in the plane, which will cost an energy of order 
J ~J± per singlet, as long as spin excitations are gapped. 

Now, suppose a hole (or a projected electron) hops on 
to the adjacent layer, breaking a singlet, which creates 
two unpaired spinous, one in each plane. The cost in 
energy is f2, which can not be removed by the hole 
hopping on to a third layer since that would create more 
unpaired spinous costing more energy. Now, there are 
only two processes by which the system can relax. The 
hole can hop back to the original layer, which means that 
the normal state is two dimensional. The second one is 
the two-hole process similar to the one discussed earlier, 
that is, a second hole follows the first. As in the plane, 
the pair hopping is accompanied by a singlet backflow. 
Therefore superconductivity is three-dimensional, which 
explains the 3D enhancement of Tc- 

Although this interlayer holon pair hopping is nomi- 
nally simlilar to the interlayer pair tunneling mechanism 
[59| . the physics is quite different. In the ILT theory 
electron pairs tunnel. It plays the primary role; there 
is no intra-later pair hopping. The spinon spectrum is 
gapless, the interlayer hopping matrix is diagonal in kxy^ 
i.e., long-ranged. In our case, the hopping is localized 
in real space. The primary mechanism for supercon- 
ductivity is the intra-layer holon pair hopping, which 
contibutes the the main part of the condensation energy. 
The interlayer hopping process makes superconductivity 
three-dimensional and enhances T^. 

Superconductivity in the 3D system 

The second process described above contributes the 
following interplane pair hopping term 



Hi, 



-i-V[Ft.(z)F,,(z + l) + /i.c], 



(30) 



where t^s > is an effective pair hopping amplitude 
which is (<< ts), the intraplane pair hopping amplitude. 
In this paper we will treat it as free parameter. Note that 



the vector fields =< hjhi > lie on the plane. Fur- 



thermore, this term does not modify the normal holon 
spectrum at the Hartree-Fock level. 

The MF problem for the superconducting state is again 
reduced to a 2D problem with an additional pairing term. 
The equation for the order parameter (with rf-wavc sym- 
metry) now becomes 



, tanh(/3£:(k)/2) 



(31) 



where, as before, 

i?(k) = [C2(k) + |A,(k)p]i/2^ 



but with a modified gap function 



sin ky ) . 



(32) 



The equation for x is unaffected. As shown in Fig. (7) 
_Fo increases with tzs, for given values of t^, the effect 
increasing slowly with increasing x. Note that the sign 
of tzs is determined by the symmetry of the underlying 
RVB state, and is of considerable importance since, as 
seen from the figure, Fq actually decreases if the sign is 
reversed by hand (bottom line). 
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FIG. 7: 3D enhancement of superconducting condensation en- 
ergy: Order parameter vs interplane pair hopping amplitude 
tz, which is the same as t^s in the text. It is enhanced in the 
actual case (dotted line). However, it is reduced if the sign of 
the interaction is reversed by hand (bottom line). 

IX. Electron doping vs hole doping 

High- Tc superconductivity also occurs in electron- 
doped systems. This is expected theoretically since the 
square-lattice Hubbard (or t-J) model has an electron- 
hole symmetry about half filling. The symmetry is re- 
alized by the transformation: Cia — > c^^ on one sublat- 
tice and Cia — *■ —c]^ on the other. If we assume that 
the orbital structure in the Cu02 plane is roughly the 
same for the two cases, an electron-doped system at an 
electron concentration 1 + x should exhibit roughly the 
same behavior as a hole-doped system at a concentration 
1 — x. Experimentally this is not found to be the case. 
In the electron-doped system the superconducting state 
occupies a much smaller region in the parameter space, 
and Tc's are also much smaller. The AF insulator phase 
occupies a much larger region. 

In a real system there will always be an asymmetry 
since one has to include an intra-sublattice hopping term 
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in the original model: 

H' = t'Y,clcia + h.c., (33) 

il 

where i and I are next-nearest neighbors, and this term 
changes sign under the e-h transformation. A rough esti- 
mate from ARPES and electronic structure calculations 
is: t'/t « (0.1 - 0.3) [63|. The t' term would generate 
a J' = 4t'^/J7, which is very small since J' / J = {t'/t)'^, 
and can be neglected. In general, t' should not change 
the physics qualitatively; and for a Fermi-liquid or a nn 
state the quantitative difference should not be too large 
since physics is dominated by the t term. The observed 
difference is much larger than what is expected. 

In our case the difference can be significant since t is 
renormalized away, and the hopping parameter th that 
characterizes renormalized intra-sublattice hopping (see 
Eqs. 17-18) is only ~ J. It should be noted this term does 
not change under the e-h transformation. Therefore, the 
effect of the t' is to increase the net intra-sublattice hop- 
ping amplitude in one case and decrease it in the other. 
Of course, t' will also be renormalized to t'^ff, but RG 
effect will be smaller since this term moves an electron 
within the same sublattice which has ferromagnetic cor- 
relations. Now, in the MF approximation H' becomes 

^' = -t'effT.^'2B.jh]h. + A, E ^L^- + '^■^■l' (34) 

ij a 

which generates intra-sublattice holon hopping. Here 
Bij =< b\^bja > and Dij =< h^jhi > are the average 
hopping amplitudes for spinons and holons along the di- 
agonal. Above T* this term does not contribute since 
Bij = = Dij . We assume that the t'^j^j: is small enough 
that it does not by itself break gauge symmetry and gen- 
erate coherent hopping of holons above T* . 

The situation is quite different below T*, since the di- 
agonal holon hopping amplituc Dij ^ in the pseudo- 
gap metal; and the diagonal Bij is also nonzero in the 
presence of the singlet condensate. The symmetry of 
B*j is easily calculated from the MF theory, and we find 
B*j = along the diagonal (1, ±1) directions, where 

B2 = smk,smky^[l/2 + niujik))]. (35) 

Since the spinon spectrum has minima at k = (7r/2, 7r/2), 
B2 > 0. Thus the t' term contributes 

ejj(k) = — 4t5j sin fca, sin fcj,, (36) 

to the holon spectrum, where t'j^ = t'^f^B2- Therefore, 
the extra holon hopping term in the diagonal direction 
has the same symmetry as the earlier one; except that 
the sign is opposite for the hole-doped case for which, 
tjj > 0. As shown in Fig. 8 (upper panel), the extra 
term has maxima at ±(7r/2, — 7r/2), where the original 
spectrum is a minimum, making the band shallower and 



holons heavier (and increases the holon DOS). Since t' 
term does not contribute to the pairing interaction, the 
net effect is to reduce pair breaking, which increases Fq 
and Tc- This is seen in Fig. 9 where we plot i^o vs x at 
T = 0. 




FIG. 8: The origin of electron-hole asymmetry: the holon 
spectrum in the second quadrant in the presence of next- 
nearest-neightbor hopping t'^jj — ±0.2te//. For positive t'^ff 
(hole-doped case) the minimum becomes shallower (upper 
panel), which would to lower Fermi velocity and higher Fq. 
For the electron-doped case, t'^jj is negative, and minimum 
deepens (lower panel), which would lead to higher fermi ve- 
locity, and lower _Fo- 

The problem of an electron-doped system at an elec- 
tron density of 1 + x is the same as that of a hole-doped 
system at a hole density of x, except that the sign of t', 
(and hence, that of i'^) is now changed. In this case, the 
total holon hopping amplitude along the diagonal add. 
The extra term has minima at the minima of the orig- 
inal spectrum (Fig. 8 - lower panel), making the band 
steeper, which making holons lighter and thus decreases 
the DOS. As shown in Fig. 9, Fq is now smaller, as one 
would expect. We stress that these effects are amplified 
because the pairing interaction and holon hopping are of 
the same order of magnitude. 

Interestingly the theory also provides a straightfor- 
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FIG. 9: Electron-hole asymmetry: the condensation energy 
vs the effective diagonal hopping amplitude t'f^ for fixed x. 
The hole doped case (dotted line) corresponds to t'f^ > 0, and 
the electron doped case (broken line) corresponds to tj, < 

ward explanation for why the size of the AF region 
is larger for the electron-doped case simply on the 
basis of the symmetry of the RVB state. The point is 
that average spinon hopping amplitude Bij measures 
ferromagnetic correlations between sites i and j. It is 
nonzero when i and j are on the same sublattice, and 
has the largest value when the system has long-range 
AF order. This can be verified from Eq. (35) which 
shows that right-hand side is largest when the spinon 
energies w(k) are smallest, i.e., when the spectrum is 
gapless (corresponding to long-range AF order). We 
also see from Eqs. (20,21) that Dij = ±Di, with Di > 
along the (1,±1) directions. Hence, the contribution of 
the t' term to the energy per bond is t'^j^B2Di which 
has the same sign as t' . It follows that AF ordering will 
be opposed by this term for the hole-doped case [t' > 0), 
and will be favored for the electron-doped case (t' < 0). 

X. Conclusion: 

In this paper we have shown that assumptions of con- 
tinuity from half filling and renormalization gives rise to 
the number and type of phases that are seen in the exper- 
imental phase diagram. The theory makes three major 



predictions: a spin gap which plays a pivotal role in all 
three phases; a spinless Fermi liquid below T*, but one 
without without a Fermi surface; and a state above T* 
which is not a metal but a novel quantum insulator. Ex- 
cept for superconductivity, the motion of holons do not 
introduce any new order, other than those at half filling. 
The normal states are predicted to be rather different 
when compared with other theories. We have presented 
a careful review of the old and new experimental results 
in the charge, spin and electron sectors, and found very 
strong evidence in support of these predictions. 

Superconductivity also occurs naturally, via pair hop- 
ping of holons, driven by the existence of pre-paired 
spinous. It is found that the possible symmetries of the 
superconducting order parameter are severely restricted 
because the symmetry of the spin part of the electron pair 
is already determined at half filling. We have shown that 
this leads to a robust c?-wave. Remarkably, the theory 
provides a natural explanation for the two-dimensionality 
of the normal state, and also of the pronounced differ- 
ence between hole-doped and electron-doped cuprates. 
Most of results follow from the type and symmetry of 
the renormalized Hamiltonians, detailed calculations are 
not necessary. 

The theory provides a basis for future calculations of 
finite-temperature properties, which would require care- 
ful treatment of fluctuations, including gauge fluctua- 
tions, separately for each phase. In particular, because of 
the low dimensionality, and since the energy scales of all 
the terms in the renormalized Hamiltonian scale with J 
and are comparable, the mean-field approximation would 
not work near the transitions. For this reason, probing 
the nature of the transition near T* , and the calculation 
of Tc{x) or, possibility of a Nernst effect, would require a 
more sophisticated treatment. We also do not discuss the 
issue of collective modes, the treatment of which would 
require going beyond the MF approximations. There are 
a number of such collective excitations. An important 
one is the physical electron itself, which can be probed 
by the photoemission spectrum. Experimentally one sees 
something like a Fermi surface above T*, and a Fermi 
arc in the pscudogap normal state, although there arc no 
quasiparticle excitations associated with these. One im- 
portant issue is the possible appearance of sharp peaks in 
the gap region of the superconducting state. Our theory 
does not rule out the emergence of a more conventional 
Fermi liquid state at large x, where the physical electron 
would appear as a collective excitation via spin-charge 
recombination [2^ . These and other issues would be dis- 
cussed in a future paper. 



[1] J. G. Bednorz and K. A. Muller, Z. Phys., 64, 189 (1986). (1988). 

[2] T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999). [5] D. P. Arovas, and A. Auerbach, Phys. Rev. B 38, 316 

[3] P. W. Anderson, Science 235, 1196 (1987). (1988). 

[4] G. Baskaran and P. W. Anderson, Phys. Rev. B 37, 580 [6] S. Sarker, C. Jayaprakash, H. R. Krishnamurthy, and M. 



18 



Wen, 



1694 



Ma, Phys. Rev. B 40, 5028 (1989). 
[7] C. Jayaprakash, H. R. Krishnamurthy, and S. K. Sarker, 

Phys. Rev. B 40, 2610 (1989); D. Yoshioka, J. Phys. Soc. 

Jpn. 58, 1516 (1989). 
[8] Z. Y. Weng, D. N. Sheng, and C. S. Ting, Phys. Rev. 

Lett. 80, 5401 (1998). These authors connect with half 

filhng by using a representation in which both holons and 

spinous are bosons. But they have not shown their states 

are stable relative to spiral states. 
[9] G. Baskaran, Z. Zou, P. W. Anderson, Solid State Com- 

mun. 63, 973 (1987). 
[10] G. Kotliar and J. Liu, Phys. Rev. B 38, 5142 (1988). 
[11] L Ichinose, T. Matsui and M. Onoda, Phys. Rev. B 64, 

104516 (2001). 

[12] P. A. Lee, N. Nagaosa and X.-G. Wen, Rev. Mod. 

Phys.78, 17 (2006). 
[13] F. M. Hu, S. K. Sarker and C. Jayaprakash, Phys. Rev. 

B 50, 17901 (1994). 
[14] P. B. Weigmann, Phys. Rev. Lett. 60, 821 (1988). 
[15] P. A. Lee, Phys. Rev. Lett. 63, 680 (1989); X. -G. 

Phys. Rev. B 39, 7223 (1989). 
[16] S. K. Sarker, Phys. Rev. B 77, 052505 (2008). 
[17] N. Read and S. Sachdev, Phys. Rev. Lett. 62, 

(1989); Phys. Rev. B 42, 4568 (1990). 
[18] S. A. Kivelson, D. S. Rokhsar and J. P. Sethna, Phys. 

Rev. B 35, 8865 (1987). 
[19] D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 

2376 (1988). 
[20] S. K. Sarker, Phys. Rev. B 46, 8617 (1992). 
[21] S. K. Sarker, Phys. Rev. B 61, 8663(2000). 
[22] S. Elitzur, Phys. Rev. D 12, 3978 (1975). 
[23] L. B. loffe and A. I. Larkin, Phys. Rev. B 39, 
[24] C. L. Kane, P. A. Lee and N. Read, Phys. Rev. B 39, 

6880 (1989). 

[25] S. Schmitt-Rink, C. M. Varma and A. E. Ruckenstein, 

Phys. Rev. Lett. 60, 2793 (1988). 
[26] D. Poilblanc, T. Ziman, H. J. Schulz and E. Dagotto, 

Phys. Rev. B 47 14267 (1993); for a review see, E. 

Dagotto, Rev. Mod. Phys. 66 763 (1994). 
[27] S. K. Sarker, Phys. Rev. B 47, 2940 (1993). 
[28] F. D. M. Haldane and H. Levine, Phys. Rev. B 40, 7340 

(1989). 

[29] N. Read and B. Chakraborty, Phys. Rev. B 40 7133 
(1989). 

[30] S. K. Sarker, Phys. Rev. B 68, 085105 (2003). 

[31] D. Poilblanc and T. M. Rice, Phys. Rev. B 39, 9749; 
H. J. Schulz, Phys. Rev. Lett. 64, 1445 (1990); J. Phys. 
(Paris)50, 2833 (1989); J. Zaanen and O. Gunnarson, 
Phys. Rev. B 40, 7391 (1989; M. Kato, K. Machida, H. 



(1989). 



[32' 
[33; 

[34' 
[35; 

[36 

[37 
[38; 

[39 

[40 

[41 
[42; 

[43' 
[44; 

[45' 
[46; 

[47 



[49 
[50 
[51 
[52; 



[53' 
[54' 
[55' 
[56; 

[57 

[58; 

[59 
[60 



Nakanishi and M. Fujita, J. Phys. Soc. Jap. 89, 1047 
(1990). 

P. B. Visscher, Phys. Rev. B 10, 943 (1974). 

S. K. Sarker, C. Jayaprakash, and H. R. Krishnamurthy, 

Physica C228, 309 (1994). 

W. W. Warren et ai, Phys. Rev. Lett. 62, 1193 (1989). 
R. E. Walstedt, R. F. BeU and D. B. Mitzi, Phys. Rev. 
B 44, 7760 (1991). 

M. Bankay, M. Mali, J. Roos and D. Brinkmann, Phys. 

Rev. B 50, 6416 (1994). 

R. Stern et ai, Phys. Rev. B 50. 

M. H. Julien et ai, Phys. Rev. Lett. 76, 4238; Bobroff et 
ai, Phys. Rev. Lett. 78, 3757 (1997). 
N. J. Curro, T. lami. C. P. Slichter and B. Dabrowski, 
Phys. Rev. B 56, 877 (1997). 

H. Zimmermann et ai, J. Less-Common Met. 164-165, 
138 (1990). 

Ch. Renner et ai, Phys. Rev. Lett. 80, 149 (1998). 

N. Miyakawa et ai, Phys. Rev. Lett. 83, 1018 (1999): N. 

Miyakawa et ai, Phys. Rev. Lett. 80, 157 (1998). 

H. Takagi et ai, Phys. Rev. Lett. 69, 2975 (1992). 

K. Segawa and Y. Ando, Phys. Rev. Lett. 86, 4907 

(2001) . 

K. Takenaka et al, Phys. Rev. B 68, 134501 (2003). 

A. F. Santander-Syro et ai, Phys. Rev. Lett. 88, 097005 

(2002) . 

B. Bucher et al, Phys. Rev. Lett. 70, 2012 (1993). 

T. Ito, K. Takenaka, and S. Uchida, Phys. Rev. Lett. 70, 
3995 (1993). 

Batlogg et al, Physica C 235-240, 130 (1994). 
Y. Ando et al, Phys. Rev. Lett. 92, 197001 (2004). 
W. J. Padilla et al, Phys. Rev. B 72, 060511(R) (2005). 
J. G. Orenstein et al, Phys. Rev. B 42, 6342 (1990); S. 
Uchida et al., Phys. Rev. B 43, 7942 (1991); S. L. Cooper 
et ai, Phys. Rev. B 47, 8233 (1993). 

A. Ino et al, Phys. Rev. B 65, 094504 (2002). 

T. Ypshida et ai, Phys. Rev. Lett. 91, 027001 (2003). 
Y. Wang et al, Phys. Rev. Lett. 95, 247002 (2005). 

G. Deutscher, A. F. Santander-Syro, and N. Bontemps, 
Phys. Rev. B 72, 092504 (2005). 

H. J. A. Molegraaf et al. Science 295, 2239 (2002); A. 

B. Kuzmenko et ai, Phys. Rev. B 72, 144503 (2005). 

D. G. Clarke, S. P. Strong and P. W. Anderson, Phys. 
Rev. Lett. 72, 3218 (1994). 

P. W. Anderson, S. Chakravarty, S. Strong and A. Sudbo, 
Science 261, 337 (1993). 

E. Pavarini et ai, Phys. Rev. Lett. 87, 047003 (2001). 



